Spontaneous quantum Hall effect via thermally induced quadratic Fermi point 
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Gapless electronic systems containing topologically nontrivial Fermi points are sources of various 
topological insulators. Whereas most of these special band-crossing points are built in the electronic 
structure of the non-interacting lattice models, we show that a quadratic Fermi point characterized 
by a non-zero winding number emerges with a collinear triple- Q spin- density- wave state that arises 
from a perfectly nested but topologically trivial Fermi surface. We obtain a universal low-energy 
Hamiltonian for the quadratic Fermi point and show that such collinear orderings are unstable 
against the onset of scalar spin chirality that opens a gap and induces a spontaneous quantum Hall 
insulator as the temperature tends to zero. 



Fermi points with nontrivial Berry connection have re- 
cently been the focus of intense theoretical and experi- 
mental effort. These singular points are the momentum- 
space counterpart of topological defects in ordered me- 
dia; both are characterized by a topological invariant [ij . 
Topological Fermi points are usually robust against small 
perturbations which preserve certain discrete symmetries 
of the system. On the other hand, when a gap opens 
at a singular Fermi point, its topological nature can be 
transferred to the resultant insulating state, giving rise 
to intriguing phases such as spontaneous quantum Hall 
or quantum spin-Hall states [2|. 

The most famous Fermi point is the Dirac point. 
The Fermi "surface" of a half- filled honeycomb lattice, 
like graphene, consists of a pair of these points. Each 
Dirac point is similar to a vortex in the Brillouin zone 
(BZ) and is characterized by a topological winding num- 
ber n = ±1 [Ij. Haldane showed that the inclusion of 
complex hoppings in a honeycomb lattice opens up a gap 
at both Dirac points, giving rise to a quantum Hall state 
in absence of external magnetic field [3| . Similar ideas for 
realizing quantum Hall insulators through gapped Dirac 
points have been explored in double-exchange models on 
square, kagome, and checkerboard lattices [j4S|- By cre- 
ating and manipulating Dirac points with the aid of arti- 
ficial gauge fields, it is also possible to control topological 
phase transitions in optical lattices 0, [![ . 

Quantum Hall effect can also be induced by open- 
ing a gap at a topological Fermi point with a quadratic 
band dispersion [9[. These points are vortices carry- 
ing m ultip le topological charges n = ±2 in momentum 
space [lo[|Tl|. Because of its higher winding numbers, 
a quadratic Fermi point can decay into several elemen- 
tary Dirac points while preserving the total topological 
charge [lo|- Unlike Dirac points where interaction ef- 
fects are suppressed due to a vanishing density of states 
(DOS), quadratic Fermi points are unstable against ar- 
bitrarily weak short-range interactions due to their finite 
DOS [sl, [ill. These Fermi points are usually protected 
in lattice models by C4 or Cq point group symmetries. 
Explicit examples are tight-binding models on checker- 



FIG. 1: (Color online) (a) Noncoplanar spin order with a 
quadrupled unit cell on a triangular lattice. The magnetic 
moments at the four sublattices point toward corners of a 
regular tetrahedron (see the inset), (b) Collinear spin order 
on the triangular lattice. In the enlarged unit cell, one site 
has large spin moment 3A while the other three sites have 
small moment —A. 



board and kagome lattices, respectively [ol, [l3|. Topo- 
logical quadratic Fermi points also appear in physical 
systems such as bilayer fflaphen e Il4l4l9| , photonic crys- 
tals [20I , oxide heterostructures (2l|7and surface state of 
topological insulators [22 1. 

present the first example of 



In this paper we present tJie tirst example ot a 
quadratic Fermi point with non-zero winding numbers 
that emerges from a topologically trivial Fermi surface 
through electron interactions. This special Fermi point 
results from a collinear triple- Q magnetic ordering that is 
favoured by perfect nesting of the original Fermi surface. 
Moreover, by analyzing the stability of this emergent 
Fermi point, we demonstrate that a non-coplanar chiral 
spin order is always more stable than any other triple- Q 
ordering (collinear or coplanar) as temperature tends to 
zero. The simple reason is that the non-coplanar state 
lowers its energy by opening a gap at the emergent Fermi 
point that also leads to spontaneous quantum Hall effect. 
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To illustrate these phenomena we consider Hubbard-like 
models on 2D lattices with a Cq symmetry such as trian- 
gular or honeycomb lattices. It has been shown recently 
that the spin- density- wave (SDW) ground state of such 
models has a noncoplanar chiral order when the chemi- 
cal potential reaches the saddle-point of the band struc- 
tures (23l-[27|. This magnetic ordering has a quadrupled 
unit cell with local spins pointing toward corners of a 
regular tetrahedron; an example of such ordering on the 
triangular lattice is shown in Fig. [TJa). 

The spontaneous quantum Hall effect of the SDW state 
arises from a non-zero scalar spin chirality (S^ • Sj x S^) 
of the underlying magnetic structure [23|, |28|. When 
conduction electrons propagate through such a spin tex- 
ture, their wavefunction acquires a Berry phase which is 
equal to half the solid angle subtended by local moments 
of each elementary plaquette. The Berry flux which is 
indistinguishable from a real magnetic flux induces the 
quantum Hall effect. Indeed, quantized Hall conductiv- 
ity cTxy = /h has been explicitly calculated for the chi- 
ral SDW state in both triangular and honeycomb lat- 



tices 



Here we show that the recently discovered finite- 
temperature SDW phase with collinear magnetic mo- 
ments [Isl provides a useful starting point for under- 
standing the origin of the spontaneous quantum Hall ef- 
fect that necessarily appears at lower temperatures. A 
striking feature of this collinear SDW state shown in 
Fig. mb) is the existence of a quadratic Fermi point 
in its mean- field band structure [29]. Gapless charge 
excitations in this SDW state exist only for one spin 
branch. More importantly, we show that the quadratic 
Fermi point has a topological winding number n = ±2. 
We derive a universal low-energy Hamiltonian for exci- 
tations around the Fermi point. Remarkably, the gap 
of the quantum Hall state is proportional to the scalar 
spin chirality of the ordered spin state. Consequently, the 
transition from collinear to chiral SDW order coincides 
with the onset of quantum Hall effect. 

We start by considering SDW instabilities in 2D lat- 
tices with a Cq symmetry, such as triangular, honeycomb, 
kagome and their decorated variants [30[. The tight- 
binding DOS includes one or two Van Hove singularities. 
For example, a logarithmically divergent DOS appears 
for filling factor 3/4 in the triangular lattice, and 3/8 or 
5/8 in the honeycomb lattice. At these filling factors, the 
Fermi surface is a regular hexagon inscribed within the 
BZ; see Fig. [2fa). Remarkably, pairs of parallel edges of 
this Fermi surface are perfectly nested by wavevectors 
which are equal to half of reciprocal lattice vectors. The 
perfect Fermi surface nesting is quite robust and is bro- 
ken only by the inclusion of third nearest-neighbor or 
longer range hopping amplitudes. 

SDW order appears as a weak-coupling instability in- 
duced by perfect Fermi surface nesting. In particular, 
the nature of the SDW instability is mainly controlled by 
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FIG. 2: (Color online) (a) The Brillouin zone of lattices with 
a Cq symmetry and the Fermi surface when electron filling 
reaches the Van Hove singularity of the band structure. The 
three nesting wavevectors are Qi = (27r/y3,0) and Q2,3 = 
(-7r/V3,dz7r). (b) Phase diagram of the triangular-lattice 
Hubbard model. SDWi and SDW2 denote spin-density-wave 
states with collinear and noncoplanar magnetic moments, re- 
spectively, (c) Configurations of the three vector order pa- 
rameters A^n in the collinear (top) and chiral 'tetrahedral' 
(bottom) SDW phase. 



the states that are close to the saddle points shown 
in Fig. [2fa) {r] = 1,2,3), where a vanishing Fermi ve- 
locity gives rise to a logarithmically divergent DOS. The 
effective Hamiltonian for the SDW ordering expressed in 
terms of these low-energy electrons is [29| 



(1) 
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^/C;a;/3k;q;p 

+ X] X] ^a,/34a,k-q4/3,k+q^C/3,k+p^^a,k-p ' 
^/C;a;/3k;q;p 

where c]^^ ^ creates an electron with spin a = t, | and 
momentum K^^ + k (|k| <C 1), while g and denote 
the forward and umklapp scatterings, respectively (ioj . 
The dispersion around the three saddle points is given 
by (k) = /c^/c^/m*, where 1/m* ~ 0{t) is of the order 
of bandwidth, (r^^C) is a cyclic permutation of (123), and 

= k- e^, with 63 (1,0) and ei,2 (-|, For 
the special case of Hubbard model, we have g = g' = 
US^^^/N^ where U is the on-site Coulomb repulsion, P 
denotes the opposite spin, and N is the total number 
of sites. Although renormalization-group (RG) analysis 
indicated that superconductivity instability is asymptot- 
ically dominant 3lj, the SDW vertex is the largest at 
intermediate RG scale and becomes dominant by slightly 
doping away from the Van Hove singularity [32]. 

The effective filling fraction restricted to this low- 
energy model can be obtained by a simple geometrical 
consideration. As shown in Fig. [2ja), the filled region in 
the vicinity of a saddle point is bounded by two straight 
lines with an angle = 27r/3 between them. This gives 
rise to an effective filling factor: z/gfr = O/tt = 2/3. In the 
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folded BZ of the triple-Q SDW states, the three saddle 
points are shifted to the zone center; each corresponds 
to two electron bands with different spin species. An ef- 
fective 2/3 filling indicates that 4 out of the 6 bands are 
filled in the insulating state. 

The component of the SDW order parameter as- 
sociated with the nesting wavevector is = 

I Ek <4«,k^a/3Cc/3,k)' where v = g ^ g' . A mean-field 
decoupHng of the interaction terms in ^ based on this 
SDW order parameter yields a Hamiltonian: 



with the interaction matrix: 




As ■ a A2 • cr 

kskil Ai - a 
Ai - a kik2l 



(3) 



Here I is a 2 x 2 identity matrix, a = (<j^, a^, <j^) is a 
vector of Pauli matrices, and we have set m* = 1 for 
simplicity. 

To examine the structure of the SDW state near the 
ordering temperature, a Ginzburg-Landau expansion up 
to sixth order in A was derived in Ref. [29| by integrat- 
ing out the fermions. The analysis found that a collinear 
SDW state with Ai = A2 = A3 is favored immediately 
below the magnetic transition; the corresponding spin or- 
der is shown in Fig. [IJb). At a lower temperature, the 
system undergoes another transition into a chiral SDW 
state whose three components of the order parameter are 
orthogonal to each other and have the same amplitude, 
as shown in Fig. [2fc). The real-space magnetic order 
shown in Fig. [T]^a) has a quadrupled unit cell with lo- 
cal spins pointing toward different corners of a regular 
tetrahedron. An explicit calculation for the triangular- 
lattice Hubbard model gives a phase diagram shown in 
Fig.[2fb), which is consistent with the two-stage ordering 
scenario described above. 

The collinear SDW state is a half-metal which hosts 
a Fermi point at the center of the folded BZ; its mean- 
field band structure is shown in Fig. [3fa). To under- 
stand the nature of this gapless point, we first focus on 
the fermion spectrum at k = 0. We assume that the 
three order parameters have the same amplitude A, as 
such states are favoured by the fourth order contributions 
to the energy expansion in powers of A^ |29i]. The in- 
teraction matrix (j3j) at the F point in momentum space 
is equivalent to a tight-binding problem on a triangu- 
lar loop (see Fig. [3jc)) with a spin-dependent hopping 
Uj = t Uij = A^ • a = Aexp(z|- J A • dr), where A is a 
non-Abelian gauge potential. Its spectrum depends only 
on the gauge- invariant Wilson loop W = ti {U12 U23 Usi). 
Interestingly, this non-Abelian fiux is proportional to the 
scalar spin chirality 
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FIG. 3: (Color online) The mean-field band structure along 
k = (/c, 0, 0) in the folded BZ for the (a) collinear and (b) 
noncoplanar SDW states. The calculation is done for the 
special case of triangular lattice. The quadratic Fermi point 
with winding number n — ±2 is marked by in panel (a). 
Each band in the noncolanar SDW state is doubly degenerate, 
(c) An equivalent tight-binding model on the triangular loop 
for Hamiltonian (|3]). (d) Energy levels at the k = as a 
function of the invariant Wilson loop W . The collinear and 
tetrahedral SDW states have = and 1, respectively. 



Explicitly, the energy levels at the F point are solutions 
of the polynomial equation: e^{e^ — 3)^ +4(|W| — 1) = 0, 
and always appear in pairs. Fig. [3fd) shows the 

three positive eigen-energies as a function of the non- 
Abelian fiux. Noting that there are three bands with 
negative energy, the charge gap at 2/3 filling is given by 
Ae = £2 — ^1- This gap closes only when VF = 0, i.e. for 
any coplanar, and particularly collinear, SDW state. 

We now show that the transition from collinear to 
chiral SDW states is a topological phase transition in- 
volving a quadratic Fermi point. Without loss of gen- 
erality, we assume Ai = A2 = A3 = Az for the 
collinear order. The low-energy electrons in the vicinity 
of the emergent Fermi point come from the two bands 
e\ and £2 in Fig-ISl^d); the corresponding eigenstates are 
V^i = ^(cu-C2;) and V^2 = ^(ci; + C2; - 2c3;), respec- 
tively. Introducing a pseudospin = ±1 to label these 
two states, the low-energy Hamiltonian in this basis is 
given by 



(5) 



W 



-iAi • (A2 X As)/ A' 



where and are Pauli matrices and 
ho = A-ikl + kl)/4, 
= d, = (kl - kl)/A. (6) 

A diagonalization of 'H{k) yields a flat band ei^k = A and 



(4) a quadratic band e2,k = A — |kp/2. The pseudovector 
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field d = {dx^ d^) shown in Fig. HJ^b) lias (i-wave symme- 
try. The topological charge for the singular k = point 
is given by the winding number of pseudovector field: 
n = ^ V6>(i(k)-(ik = ±2 where 6^ = diTctdin{dz/dx) 
and C is a contour enclosing the Fermi point. 

The existence of this topological Fermi point is pro- 
tected by the Cq symmetry preserved in the collinear 
SDW phase and the spin collinearity masqueraded as 
an effective time-reversal symmetry. Further symmetry- 
breaking at lower temperatures could remove the 
quadratic Fermi point by either splitting it into elemen- 
tary Dirac points or simply opening a gap. To study the 
stability of the collinear SDW state, we introduce small 
deviations to the order parameters: = Az + m^ with 
± z. Substituting into Eq. ([2]) and projecting to the 
low-energy doublet manifold ^ = (V^i, ^^2), the mean- field 
Hamiltonian becomes 



.\/6A 



\/3A2 



(7) 



*(r), 



where 1/v' = 1/'U + 1/6A is the effective inverse coupling. 
The doublet order parameter (Qi, Q2) given by 



Qi = (|mi|V|m2|'-2|m3|')/V6, 



Q-2 



|mi| 



Imal 



(8) 



describes a nematic phase for the ^ fermions in which 
the Cq rotational symmetry is broken down to C2 by 
splitting the quadratic Fermi point into two Dirac points. 
The corresponding SDW order is dominated by a single 
nesting wave vector. The nematic phase remains a half- 
metal. Since the total winding number is conserved, the 
two residual Dirac points in the nematic phase carry the 
same topological charges. The order parameter n is the 
scalar spin chirality (j4]): 



Ai • (A2 X A3) = A(mf ml - m\ (9) 

+ 7712 ^3 ~ ^2 ^3 + ^3 ^1 ~ ^3 ^1 ) ' 



It characterizes an insulating phase with broken time- 
reversal symmetry and a zero-field quantized Hall con- 



ductivity (7x 



Minimization of Eq. (O yields a ground state that pre- 
serves the Ce rotational symmetry while breaking the ef- 
fective time-reversal symmetry. The transition from the 
/i: = phase into the quantum Hall state with 7^ is a 
discontinuous one. To see this, we compute the ground- 
state energy of ([7]) as a function of k [33|: 



V 



6A4 k|' 



(10) 



where p* is the density of states at the quadratic Fermi 
point, A is an ultraviolet cutoff, and V is the system 
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FIG. 4: (Color online) (a) Fermi point at a quadratic band 
crossing at k = 0. (b) The pseudovector field d = (dx^dz) 
with dx — (kx — ky)/4 and dz — kxky/2 resembles a vortex in 
XY systems with a winding number n = 2. 



volume. The function E{i<i) has two minima si hi = {) 
and K ~ A/y^. At higher temperatures, the system 
may remain in the = minimum corresponding to 
the collinear SDW state, as it was found in Ref. [29] for a 
Hubbard model on triangular and honeycomb lattices. 
Upon lowering the temperature, the system switches 
to the absolute minimum with ti = ±A^ (since the 
triple product is bounded) which corresponds to the non- 
coplanar tetrahedral order. 

The first-order transition scenario is consistent with 
our explicit mean-field calculations of the Hubbard model 
on triangular or honeycomb lattices. It is worth noting 
that the collinear SDW phase is stabilized by thermal 
fluctuations at finite temperatures. A lower free energy 
of the uniaxial SDW state results from the gapless elec- 
tronic excitations and the softer magnetic fluctuations as- 
sociated with collinear spin order. Although the collinear 
SDW state hosting a topological Fermi point has been 
shown to exist in three of the most representative 2D lat- 
tices with Cq rotational symmetry [29|, [33| , and maybe 
also in their decorated variants, its stability at finite tem- 
peratures depends on details of the model. 

To summarize, we have shown the emergence of a 
quadratic Fermi point characterized by a winding num- 
ber n = ±2 in triple-Q collinear SDW states of different 
two-dimensional lattices. Such SDW phase arises from 
a perfectly nested Fermi surface in lattices with a six- 
fold rotational symmetry such as triangular, honeycomb, 
or kagome. Unlike most topological Fermi points which 
result from the special band structures of certain lat- 
tice problems, this quadratic Fermi point emerges from 
a topologically trivial Fermi surface through electron in- 
teractions. We have also obtained a universal low-energy 
Hamiltonian for the quadratic Fermi point and showed 
that the instability towards a quantum Hall insulator is 
described by an order parameter which corresponds to 
the scalar spin chirality. Our theory thus explains why 
the topologically non-trivial (non-coplanar) SDW order 
is always more stable than the topologically trivial copla- 
nar states at low temperatures. 
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Although spontaneous quantum Hah effect has been 
extensively discussed in the context of double-exchange 
or Kondo-lattice models, most of these studies assume 
classical localized moments to begin with. It is then nat- 
ural to ask if this phenomenon survives for models that 
incorporate quantum spin fluctuations. Our analysis of 
the noncoplanar SDW ordering in Hubbard-like models 
provides a positive answer to this important question. 
Possible material realizations of the collinear SDW phase 
that hosts the topological Fermi point include triangular 
compound Nao.5Co02 [1^ and doped graphene (sHj . 
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SUPPLEMENTARY MATERIALS 
KAGOME LATTICE 

The triple- Q spin- density- wave (SDW) states in tri- 
angular and honeycomb lattices have been discussed in 
Refs. Here we consider another representative two- 



dimensional lattice with Cq rotational symmetry: the 
case of kagome lattice. The existence and stability of 
non-coplanar tetrahedral SDW order on kagome lattice, 
shown in Fig. Efa), has been demonstrated in Ref. 
The non-coplanar tetrahedral order shown in Fig. [5ja) 
can be described as 



Sa{r) = +Aicos(Qi -r) + A2Cos(Q2 -r) + A3Cos(Q3 -r), 

S(3{r) = +Aicos(Qi • r) - A2Cos(Q2 • r) - A3Cos(Q3 • r), (11) 
S^(r) = -Ai cos(Qi • r) - A2 cos(Q2 • r) + A3 cos(Q3 • r), 



where a, /3, and 7 denote the three sublattices, and vector 
order parameters A 1,2, 3 have the same amplitude and are 
orthogonal to each other. Since the unit cell of kagome 
lattice contains three inequivalent sites, there are several 
triple- Q collinear SDW orderings with different arrange- 
ments of spins within the unit cell. The collinear SDW 
state that is topologically connected to the tetrahedral 
order is obtained from Eq. (1) above by continuously 
closing the solid angle between the three vector order 
parameters A^. By setting Ai = A2 = A3 = Az, the 
resulting collinear SDW state is shown in Fig. [5fb) in 
which 3 sites in the quadrapuled unit cell have magnetic 
moment 3A while the moment at the remaining 9 sites 
is -A. 

We now consider the Kondo-lattice model, which 
can be viewed as the mean-field approximation of the 
Hubbard-like model, on the kagome lattice: 



H 



JH^^i-c\o,(Ta^Cj^^. (12) 



where = \{c\ ^(Tc^fBCj^f^) is the local magnetic moment. 
The band structure of the nearest-neighbor tight-binding 
model {Jh = 0) is shown in Fig. [6]^a). A logarithmically 
divergent density of states occurs at 1/4-filling due to 
perfectly nested Fermi surface. Figs. Mjo) and (c) shows 
the mean-field band structures in the presence of collinear 
and tetrahedral SDW order, respectively. At 1/4-filling, 
the lowest 6 bands are completely filled, giving rise to a 
quadratic Fermi point at k = (F point in the folded 
Brillouin zone) [Fig. [6l^b)]. On the other hand, the non- 
coplanar tetrahedral order opens a charge gap as shown 
in Fig. [6fc). 

With the addition of the kagome case discussed here, 
we have demonstrated the existence of triple- Q collinear 
SDW state and its topological connection to the non- 
coplanar SDW order in three of the most representative 
2D lattices with Cq symmetry, namely triangular, hon- 



eycomb, and kagome lattices. Although perfect Fermi 
surface nesting and divergent Van Hove singularities also 
exist in other lattices such as the decorated variants of 
the above three lattices, the properties of the triple-Q 
SDW ordering will be left for future studies. 



GROUND-STATE ENERGY 

In this section we discuss the derivation of ground state 
energy Eq. (10) in the main text. 



E{k.) 2\i^\ 



6A4 ^""^ k|' 



(13) 



The first term is the energy cost when spins deviate from 
the collinear limit. It comes from the term ^ l^^r^P 
in the effective Hamiltonian Eq. (7) of the main text. As 
also discussed there, the ground state preserves the 
rotational symmetry, hence the nematic order parame- 
ter vanishes Qi = Q2 =0. The scalar spin chirality 
is then maximized when the three vectors lie in the 
same plane perpendicular to z and point at 120° to one 
another. This gives rise to l^ryP = 2\k,\/V3A. 

The presence of spin chirality opens a gap at the 
Fermi point. The dispersion of the two electron bands 
near the F point is: 61^2 (k) ~ ho{k.) ± Y^|k|^ + J^, where 
the energy gap S = \k,\/^/3A'^. Standard calculation gives 
an energy gain J (Pk{^/\k\^ + (5^ - |kp) ~ k'^ \og{A/\K\) 
which is logarithmically divergent with tz. Here A denotes 
an ultraviolet cutoff. This accounts for the second term 
in Eq. (3). 

Fig. [3 shows the ground state energy E{i<i) as a func- 
tion of the scalar chirality for a generic case. The first 
term in Eq. (3) which is linear in |/^| gives rise to a lo- 
cal minimum at = 0, corresponding to the collinear 
SDW state. The function E{k:) has a global minimum 
8it K, ^ A/^ye coming from the second term in Eq. (3). 
However, since the scalar chirality is bounded |/^| < A^, 



FIG. 5: (Color online) (a) The non-coplanar arrangement of local moments S(r) on the kagome lattice. Different colors 
correspond to the four non-coplanar spin directions pointing toward the corners of a regular tetrahedron, (b) The collinear 
triple-Q SDW state on kagome. The local moments at the and sites are 3A and —A, respectively. The quadrupled unit 
cell in both SDW states contains 12 sites. 



(a) (b) (c) 




FIG. 6: (Color online) (a) Tight-binding band structure of kagome lattice shown with the folded Brillouin zone, (b) and (c) 
are the mean-field band structures of kagome Hubbard model in the presence of collinear and non-coplanar (tetrahedral) SDW 
ordering, respectively. At filling fraction v — 1/4, there exists a quadratic Fermi point at the F point in the collinear SDW 
state, while the spectrum is gapped in the non-coplanar SDW state. 



the ground state is reached when the three order param- 
eters are perpendicular to each other, corresponding 
to the non-coplanar tetrahedral SDW order. 
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FIG. 7: (Color online) Ground-state energy E(k) as a function 
of the spin scalar chirality k, = Ai • (A2 x A3). 



